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Summary. At present the ground-based Very High Energy (VHE) gamma-ray 
astronomy is racing to complete construction of a number of modern gamma-ray 
detectors, i.e. CANGAROO III, MAGIC, H.E.S.S., and VERITAS. They should be 
fully operational in a year's time. After much debate, the further development of 
this gamma-ray astronomy in the foreseeable future must be widely anticipated to 
proceed with the designing and building of a new instrumentation, which is primar- 
ily intended for the further drastic reduction of the energy threshold in gamma-ray 
observations down to about 10 GeV. On the ground one can hardly reach such 
low energy thresholds without considerably larger, up to 30 m diameter, optical 
telescopes, which might be able to collect sufficient amount of Cerenkov light from 
the atmospheric gamma-ray showers of that low energy. If not taken off the ground 
entirely (like GLAST), then it seems to be profitable to mount future low energy 
Cerenkov telescopes at higher altitudes in the atmosphere in order that they will 
be able to detect substantially more unabsorbed Cerenkov light from a shower. 
There are a few sites up on the high mountains of roughly 5 km height worldwide, 
which can be used for such a venture. However, one has to remember that actual 
time profile, and in particular a two-dimensional distribution (image) of Cerenkov 
light flash from atmospheric gamma-ray showers, undergo a rapid change after an 
increase in the observational level. This paper briefly describes the results of a 
topological analysis of Cerenkov light images calculated at both conventional and 
desirable altitudes of 2.2 and 5 km above the sea, respectively. A discussion on 
major upgrades of image topology at high altitude is also given. 



1 Introduction 

The usual way to proceed with design studies for a future project, at least 
in a field of VHE gamma-ray astronomy, is to perform full scale simulations 
of detector response for various anticipated event types. Apparently, gamma- 
ray-induced atmospheric showers represent a sample of signal events, whereas 
cosmic-ray showers form a dominating background (for review of Cerenkov 
imaging technique see For a future low-threshold instrument a correct 
tuning of detector design to optimize its response to gamma-ray events be- 
comes a most important issue, due to the fact that the sensitivity of such 
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a detector will be given by angular and energy resolution for signal events. 
Both of these will finally determine an efficiency of background rejection. 
Here we have studied how a single parameter of detector design, particularly 
the height of the observational site, may affect a topology of signal events. 
Therefore, we have calculated the response of a ground-based Cerenkov light 
telescope of 30 m diameter, assuming an ideal focal-plane detector. Compara- 
tive analysis of simulated events helped us to understand what are the major 
differences in parametrization of the Cerenkov light flashes from gamma-ray 
induced air showers, registered at two observational levels, i.e. 2 and 5 km 
above sea level. Considering the analysis results, there is a discussion at the 
end of this paper as to which observational level might be considered as more 
favorable for effective shower imaging. 

2 Simulations 

Numerous comparisons of a few shower simulation codes have been recently 
performed by different groups. They have revealed a rather good level of 
agreement in basic parameters of Cerenkov light emission in gamma-ray- 
induced air showers over broad energy range, starting from 100 GeV and 
expanding up to 20 TeV. Calculations presented here have been carried out 
using one of those simulation codes, namely ALTAI code 0. This numeri- 
cal code was extensively used for production of the simulated data for the 
HEGRA system of five imaging atmospheric Cerenkov telescopes at La Palma 

Shower simulations have been done for the standard continental atmo- 
sphere (US standard atmosphere) 0] for the wavelength range of Cerenkov 
light photons from 300 to 600 nm. Absorption of Cerenkov light photon in the 
atmosphere due to Rayleigh and Mie scattering was modelled according to the 
data given in gj |3] . The detector simulation procedure used here accounted 
for all efficiencies involved in the process of the Cerenkov light propagation 
and registration [3]. It includes (i) mirror reflectivity; (ii) the acceptance of 
the funnels placed in front of the photomultipliers (PMTs) (iii) the photon- 
to-photoelectron conversion inside the PMTs (bi-alkali photocathode) . 

Shower simulations were made here in the so-called "batch" mode. A 
shower propagating time, corresponding atmospheric depth, and a number of 
emitted Cerenkov photons were saved for each multiple-scattering segment 
of all electron trajectories in a shower. The actual segment size was chosen 
as small as 0.1 [gr/cm 2 ]. CPU time needed for simulation of such low en- 
ergy gamma-ray shower (E o =10 GeV), using customary computers, is short, 
and it is not an issue for any scheme's optimization. One record for a single 
multiple-scattering segment was treated as an individual "batch" of emit- 
ted Cerenkov photons. At the second step of this simulation procedure all 
recorded batches were restored and finally used in the calculation of the re- 
sponse of a number of Cerenkov telescopes, situated at different atmospheric 
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altitudes. Such an approach provides an opportunity to use exactly the same 
simulated showers for various telescope arrangements at different observa- 
tional levels. It should be noted that the estimated statistical error for the 
parameters of the Cerenkov light emission given below is < 10%. 

For the next generation of ground based Cerenkov telescopes a dish of 
roughly 30 m diameter is foreseen. Issues around the design of such a big 
reflector are addressed in [Hj. The simulation setup here consisted of 12 such 
telescopes, which were arranged along one line at distances from to 300 m 
from the shower axis. Note that the same showers were used in calculations 
for each of these telescopes. It allows for a direct comparison of the telescope's 
responses at different distances from the shower axis, and the accurate study 
of fluctuations in Cerenkov light flashes at different shower impacts. 



3 Results 

Distribution of Cerenkov light emission from atmospheric showers can be 
characterized using a smooth function 

V = V {t,r,0), (1) 

which gives a mean number of photons per unit area arriving at the obser- 
vational level at time t, with space, r = {x,y}, and angular, 9 = {0 x ,O y }, 
coordinates, calculated with respect to the shower axis. The presentation 
CJ presumes averaging over a number of photons at any local spot, because 
Monte Carlo simulations have provided the list of individual photons with 
their coordinates. In the ideal case a space-angular distribution of Cerenkov 
photons in the observational plane is simply a sum of 5-functions constructed 
for each individual photon. 

The lateral distribution of Cerenkov photons at the observational level, 
which is supposed to be perpendicular to the shower axis, is given by 

p(r)=p(\r\)= f f V (t,r,9)dtdf2. (2) 

The function p(r) is the density of Cerenkov light (the number photons hit- 
ting the unit square at r). 

In a similar way one can derive a temporal distribution of a Cerenkov light 
pulse, and a two-dimensional angular distribution (image) of a Cerenkov light 
flash 

p(t) = [ [ t](t,r,e)rdr dQ Q , (3) 

JA(r a ) J n o 

q(6)= f f v (t,r,0)rdr dt, (4) 

JA(r ) Jo 
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where A(r ) and fl a are the area of the reflector and the angular camera size, 
respectively, of the telescope placed at r — r a . 

The function n(t,r,6) can be well described by a set of functions p(r), 
p(t), and q{6), given for a number of telescope locations, r$, i = 1, n. 

3.1 Cerenkov light density 

It was emphasized in [7] , that a substantial increase in Cerenkov light density 
at high altitudes in the atmosphere might be very promising for a further 
reduction of the effective energy threshold of a telescope array, which can be 
erected at a height of about 5 km above the sea. One can see in Figure^that 
for a 10 GeV gamma-ray-induced atmospheric shower the density of Cerenkov 
light at 5 km above sea level, in a range of distances of the telescope to the 
shower core limited by r < 100 m, is about a factor of 2.5 higher than 
the corresponding density at 2 km altitude. At the same time in a range 
of relatively large impact distances, r > 125 m, the Cerenkov light density 
remains the same at both observational heights, (see Figure^. 

The atmospheric depth of the shower maximum can be estimated as 
X m ax = t ln(E /E c ), where t is the radiation length in air (t a ~ 37 g/cm 2 ), 
E a is the primary energy of air-shower, and E c is the so-called critical energy 
(E c ~ 80 MeV). Thus for a 10 GeV 7-ray- induced air-shower the atmospheric 
depth of its shower maximum is about 180 g/cm 2 . A substantial fraction of 
Cerenkov light photons emitted from the shower maximum will be absorbed 
while traveling down to the observational level. Between 30% and 16% (in 
the wavelength range of 300-600 nm) for the heights of 2.2 and 5 km above 
sea level, respectively. At the same time the Cerenkov light pool shrinks sig- 
nificantly at higher altitudes. The approximate radial size of the light pool 
at 2.2 km is about 130 m (see Figure^]), whereas at 5 km it might be roughly 
limited by ~ 90 m. It results in a corresponding increase of Cerenkov light 
density by approximately a factor of 2. This geometrical effect has a major 
contribution on the increase of Cerenkov light density at high altitude as 
shown in Figure ^ 

For a large telescope of a 30 m diameter, a requirement of a minimum 
number of 15 ph.-.e. 1 in the Cerenkov light flash, which is sufficient to trigger 
the telescope, will limit the allowed range of impact distances for 10 GeV 
gamma-ray showers to roughly r < 300 m (see Figure It corresponds 
to the same effective detection area of S = 3 x 10 5 m 2 at both observa- 
tional heights. On the other side, assuming that the Cerenkov light density 
scales with primary shower energy as p oc E , one can roughly estimate the 
minimal primary energy of the gamma-ray shower, which has still sufficient 



This requirement of a minimum number of ph.-e. for the telescope trigger is in 
fact not a generic value, and it might be slightly different in certain circumstances, 
but it does not affect the general discussion given here. 
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Fig. 1. Density of Cerenkov light, p(r) [ph. — e./m 2 ], in a 10 GeV gamma-ray- 
induced atmospheric shower at two observational levels of 2 and 5 km above sea 
level, respectively. The density was measured in photoelectrons [ph. — e./m 2 ]. A 
photon-to-photoelectron conversion efficiency of x = 0.1 was assumed here, p [ph.— 
e./m 2 ] = x ■ P [photon/m 2 ]. 

amount of light at the density profile plateau (r<125 m) and can still trig- 
ger the telescope. Simple calculations yield, accordingly, a factor of 4 and 10 
for 2 and 5 km altitudes, respectively. It means that one can catch gamma- 
ray events of energy > 2.5 GeV and > 1 GeV at 2 and 5 km observational 
height, respectively, using the same telescope. It is important to mention that 
all these extremely low-energy events will be concentrated within a radius of 
roughly 100 m, which is determined by the actual shape profile of a lateral 
distribution function of Gerenkov light. The drastic drop in photon density 
beyond 120 m will prevent the detection of such gamma-ray showers at larger 
distances to the shower axis. As a result the detection area for these events 
will be a weak function of primary shower energy. Furthermore, the detec- 
tion area of high energy gamma-ray showers will be of the same size at both 
observational heights. 

3.2 Time profile of Cerenkov light flash 

The longitudinal extension of an atmospheric shower, and its location in 
space with respect to the telescope, finally determine a time profile of the 
Cerenkov light flash. For a gamma-ray shower of a certain primary energy, 
e.g. of 10 GeV, recorded at fixed observational level (e.g. 2 or 5 km above 
sea level), the shape of the Cerenkov light flash basically depends only on 
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Table 1. Parameters of the fit in Eq.|S] D [gr/cm 2 ] denotes an atmospheric depth at 
a given observational level. R is the impact distance of the telescope to the shower 
axis. 
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the actual distance of the telescope to the shower axis. The arrival time of 
Cerenkov photons onto reflector is t — t e +tQ, where t e and is accordingly a 
propagation time of emitting electrons and photons, respectively. Electrons of 
an energy, which is sufficient for emission of Cerenkov light in the atmosphere, 
are apparently moving faster than the emitted photons. Thus at relatively 
small distances of the telescope to the shower axis (r <120 m) Cerenkov 
photons, emitted at later stages of the shower development in the atmosphere, 
are actually arriving earlier than the photons emitted at the beginning of 
shower development. It swaps over at a large distance from the shower axis, 
because of the rather long travel path in dense atmosphere for the photons 
emitted by electrons out of a dying particle cascade 2 . Therefore, in general, 
the larger the distance of the telescope to the shower axis, the broader the 
corresponding Cerenkov light pulse. 

Time pulses of Cerenkov light flashes always show a very steep rising edge, 
and a slow fall-off. They can be well fitted by 

p{t) = C-t a {l + {t/t f). (5) 

Parameters of the fits of time pulses simulated for two observational levels 
and for a number of impact distances are given in Tabled The normalization 
condition was J^° ns p(t)dt = 1. The contribution of Cerenkov photons delayed 
by longer times than 50 ns is negligible. 

Shape of time pulse can be characterized by a few parameters, e.g. 
h = ho%— *io%) h = t 50 %-t 10 %, t 3 = t go %—tio%, where t w %, t 30 %, t 50 %, and 
^90% gi ye the time tags, which are defined as, e.g. J o 10% p{t)dt = 0.1. Results 
of calculations are shown in Figure |21 One can see that for a 500 [gr/cm 2 ] 
observational level time pulses are substantially broader at impact distances 

2 This is a well-known effect, sometimes called the "sea-gull" effect, in the shape 
of Cerenkov light pulses (e.g. see [HI) 
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Fig. 2. Parameters of Cerenkov light time pulses, ti (1), t% (2), t% (3), calculated 
at two observational levels of 2 (solid curves) and 5 km (dashed curves) above sea 
level, as a function of impact distance of the telescope to the shower axis R. 

beyond 100 m. Hereafter we are dealing with Cerenkov light images averaged 
over a sample of simulated 10 GeV gamma-ray showers. For an impact dis- 
tance of ca. 250 m, the time pulse of Cerenkov light flash recorded at high 
altitude in the atmosphere will be a factor of 2 broader than for the same 
impact distance at a conventional altitude of 2 km above sea level (see Fig- 
ure 0J). Integration over the time pulse yields a total number of Cerenkov 
photons in a flash. Therefore, for a given flux of night sky background light, 
a signal-to-background ratio might be correspondingly lower by factor of 2 
for a high altitude site. For high energy gamma-ray showers (E D >100 GeV) 
Cerenkov pulses recorded at 5 km above sea level might be as long as 50 ns. 
Registration of these pulses will occur in the regime highly dominated by 
night sky background. 

It is worth noting that, at the time of writing, there is no well estab- 
lished altitude dependence of a flux of night sky background light available. 
In general this parameter is considered to be very specific for each individual 
observational site. Apparently the high altitude sites provide substantially re- 
duced attenuation and consequently more starlight from the individual stars, 
which is in fact a background for the Cerenkov telescopes. However, the effect 
of bright stars might be diminished simply by switching off the high voltage 
for those camera pixels (PMTs) collecting direct star light. Such procedure 
usually runs in automatic mode while taking the observational runs. At the 
same time one might expect a significant increase in flux of background light 
photons within the UV wavelength range (200-300 nm). However, conven- 
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Fig. 3. Time profile of Cerenkov light flashes simulated for impact distance of 
150 m and two observational levels of 2 (1) and 5 (2) km above the sea. Primary 
energy of simulated air-showers is 10 GeV. Label a.u. along Y-axis stays for the 
arbitrary units. 

tional imaging cameras are not very sensitive in a wavelength range well 
below 300 nm. 

Even though there are good reasons to believe that this flux will, in fact, 
be much lower (on average) for high altitude sites, it still needs to be measured 
at any chosen observational site. 

For the sake of thoroughness one should mention that the high altitude 
sites will noticeable increase the probability that the ionizing particles, such 
as atmospheric electrons and muons of low energy cosmic rays, can directly 
hit the camera PMTs. It will lead to some random increase in background 
light over the camera pixels. However, dedicated calculations are needed in 
order to quantify this additional component of the background, which fall 
unfortunately out of the area of this paper. 

3.3 Images 

A two-dimensional distribution of the Cerenkov light intensity in the tele- 
scope's focal plane (image), q(9), can be effectively used to derive detailed 
information about shower orientation and shape. Phenomenology of Cerenkov 
light images was discussed in [5]. For an ellipsoid-like image, the orientation 
of its major axis constrains the shower orientation in space with respect to the 
telescope optical axis. Images recorded at relatively small impact distances 
from the shower axis (R<100 m) have circular shape, and an accurate de- 
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termination of the major axis is quite difficult. For impact distances beyond 
f 00 m the image ellipsoid has a well defined shape, and the ratio of its angu- 
lar size measured along a major axis to the corresponding angular size of the 
image measured along a minor axis is ca. 2 and above. At the same time at 
large impact distances (R»200 m) the total number of Cerenkov photons in 
an image becomes rather low and high fluctuations prevent accurate measure- 
ment of the image orientation. Those two effects finally constrain the range of 
optimum impact distances for effective shower reconstruction. As mentioned 
above, the advantage of the high altitude site is mainly associated with an 
enhancement in Cerenkov light density at small distances to the shower axis 
(R<100 m). However, in this range of impact distances Cerenkov light images 
tend to have poorly determined orientation. 

Average Cerenkov light images for two observational heights of 2 and 
5 km, respectively, are shown in Figure^] Note that scales used along X- and 
Y-axis are not identical. Detailed comparison of those images revealed two 
major differences in their shape. 

1. Images of fO GeV gamma-ray showers recorded at high altitude above 
sea level are substantially larger in size (see Figure 01 and Table |2J). This 
is mainly because showers are located at a relatively smaller distance to 
the telescope, than for a conventional observational site of 2 km above 
the sea. Calculations show that the angular size of an image measured 
along the major axis, the image length, increases considerably faster with 
altitude than the angular size of image measured along a minor axis, the 
image width. This can be easily understood by comparing the ratio of 
the actual scale of a shower longitudinal development over the distance 
of shower maximum to the telescope, which is located at two altitudes 
of 2 and 5 km above the sea, respectively. In a simplified toy model the 
length of the image will scale with height of observation level above sea 
level, H , as 

L oc (H max — H D ), R « (H rnax — H Q ). (6) 

2. Images recorded at high altitude must have considerably larger displace- 
ment from the center of the focal plane. Coming closer to the shower 
maximum (H max is the height of the shower maximum) , the shower will 
be apparently seen in Cerenkov light at larger angle 

Qoctan-^ R ) (7) 

with respect to the optical axis. 

One can see in Figure ^ that the density of Cerenkov photons from a 
10 GeV gamma-ray shower is approximately the same at both observational 
heights for impact distances above 150 m from the shower axis. Thus, at 
large impact distances, the difference in image shape mentioned above is of 
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Fig. 4. Images of Cerenkov light calculated at different altitudes above sea level for 
a number of impact distances of the telescope to the shower axis as indicated in the 
pictures. Images were averaged over a sample of 10 GeV gamma-ray showers. The 
contour plots were drawn starting from the position of maximum intensity with the 
isoline increment of ln2. 



Table 2. Area A, [deg 2 ] and effective size, r = \fab [deg], (in this Table both 
values are given in a format of A/r ) of Cerenkov light images calculated for two 
observational levels of 2 and 5 km, and for a few impact distances. 
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purely geometrical origin, which is independent of the image size (the total 
number of photons in the image). At first glance, large images at high alti- 
tudes might offer better resolution for a camera of crude pixellation. However, 
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for a 10 GeV gamma-ray shower images will always be significantly affected 
by contamination of background light and reduction of background light per 
camera pixel. Large images recorded at high altitudes yield considerably lower 
Cerenkov light density per 1 str for fixed image's size. For the flux of night 
sky photons, as measured at conventional observational level of 2 km above 
the sea, the image of a 10 GeV gamma-ray shower recorded at high altitudes 
will have higher contamination of background photons per camera pixel of any 
size. 

As mentioned above, the images recorded at high altitudes must have large 
displacements from the camera center. This issue stringently constrains the 
design of the camera for the telescope placed at high altitudes, in particularly 
the field of view has to be larger. Cameras of a narrow field of view will be 
drastically limited in their ability to detect gamma-rays at high energies. 

3.4 Time-dependent imaging 

One can try to suppress contamination of night sky background light in an 
image by using a very narrow time gate. This may be tuned exactly, for 
example, to a rising edge, maximum, or tail of a time pulse. However, it leads 
to a trade off between losing a substantial fraction of Cerenkov photons and 
on the other hand a severe reduction of night sky background. This approach 
is illustrated by the images shown in Figure |SJ Photons emitted at the very 
beginning of the shower development in the atmosphere, which are mainly 
arriving at the front of the time pulse of a flash, form the so-called image 
"conk" 3 . The image of shower electrons, which propagate further into the 
atmosphere, shifts further away from the center of telescope's focal plane. 
Multiple scattering of low energy cascade electrons at later stages of the 
shower development becomes very important and it results in a very broad 
image as shown in Figure |5] One can see also in Figure that Cerenkov 
light images for such narrow time windows are in fact rather symmetric in 
shape, and they are not so good for reconstruction of image orientation. At 
the same time, an image which contains all registered Cerenkov photons 
(0 ns < t < 14 ns), has a regular ellipsoid-like shape. Note that photons, 
which are significantly delayed with respect to the front of time pulse, will 
be hitting the outer edge of the image, whereas background light photons are 
apparently dominant. 

Time-dependent imaging seems to be a rather promising approach in im- 
provement of gamma/hadron separation of extremely low energy events. De- 
spite that present analysis does not reveal an evident improvement using 
time- dependent imaging, it could be perhaps very effective in analysis based 
on centroid positions in a few triggered telescopes in an array. 



^Non-standard definition of a strongly elongated part of the comet like image, 
which has a relatively high photoelectron density. 
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Fig. 5. Average images of Cerenkov light from a 10 GeV gamma- ray shower, cal- 
culated at an observation level of 2 km above sea level, and for an impact distance 
of 150 m from the shower axis. An additional selection on arrival time of Cerenkov 
photons was applied. Corresponding time windows were 7.25-7.5 ns (a); 8.25-8.5 ns 
(b); 9-9.5 ns (c); 10-12 ns (d). Each of (a)-(d) plots contains approximately the 
same number of photoelectrons. The image shown in the lower panel was generated 
without time selection. 
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4 Conclusions 

In this paper we have attempted to perform a comparative analysis of 
Cerenkov light images simulated for two observational levels of 2 and 5 km 
above sea level, respectively. We used as the basis of our calculations a 30 m 
diameter telescope, which is a conceivable size for a future Cerenkov tele- 
scope project. The system aspect was not discussed here, but all conclusions 
are obviously relevant to any telescope of a possible future array. One has to 
mention that for a future detector, approaching a very low energy threshold 
of about 10 GeV or even below, a contamination of night sky background in 
the registered shower images is in fact a very important issue. It might finally 
constrain the choice of the observational site for such low threshold arrays of 
Cerenkov telescopes. 

Results reported here generally confirm that the observational site at 
higher altitude provides substantially higher Cerenkov light density. However 
this enhancement occurs only within the area limited by roughly 100 m from 
shower axis. Therefore all recorded gamma-ray showers well below 10 GeV 
(here we assume that the telescope has in fact sufficient area of the reflector, 
see discussion in Section 13.10 will concentrate in a region of relatively small 
impact distances. Those images are not so clearly elongated in shape, which 
makes reconstruction of the image's orientation rather difficult. 

Flashes of Cerenkov light from 10 GeV gamma-ray showers have broader 
time pulses for the impacts beyond 100 m at higher altitude. Corresponding 
images are considerably broader as well. Both of these two effects might 
apparently substantially increase the background light contamination in an 
image. Images recorded at high altitude must be further displaced from the 
center of telescope's focal plane, than for conventional altitude of ca. 2 km 
above sea level. They also, most probably, require larger field of view. 

The time-depending imaging is a very promising approach in further de- 
velopment of advanced analysis for observation of low energy gamma-ray 
showers, but it might be not very effective in resolving the problem of domi- 
nating night sky background light in the recorded images of such low energy 
gamma-rays. 

One can briefly conclude that an observational site at high altitude might 
provide further modest reduction of the energy threshold of a future detector, 
even though the shape of time pulses and in particular the topology of the 
two-dimensional angular distribution of Cerenkov light flashes recorded at 
extremely high altitudes are palpably less preferable for imaging of gamma- 
ray showers above 10 GeV. 

We hope that the results presented here may help to increase the under- 
standing of changes in topology of Cerenkov light images after an increase in 
the observational level from its conventional height of ca. 2 km up to 5 km 
above sea level. 

Ultimately, the choice of a site for the next generation of ground-based 
imaging atmospheric Cerenkov detector, which is widely believed to be a sys- 
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tern of 30 m class telescopes, will depend on both scientific and political issues 
relating to funding, international collaboration etc. The move to lower energy 
threshold is likely to remain a significant drive for the science. However, one 
has to consider all trade-offs, i.e. reasonable altitude, low level of night sky 
background, need for the robotic telescopes etc, in selecting candidate sites 
for such a detector to optimize the scientific goals. 
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